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1. Introduction 

Optical experiments at the level of single quantum emitters allow one to produce specific 
quantum states of light with photon statistics that deviate strongly from classical 
distributions [TJ |2j. Despite the experimental challenges of producing single photon 
states jHlHj, recent developments of quantum information theory have intensified interest 
in single photon sources. Realization of an efficient single-photon source (SPS) is, for 
instance, a key-problem in quantum cryptography and could more generally be applied 
to quantum information processing [5 . 

Recent experiments reported quantum key distribution (QKD) with polarisation 
encoding on single photons [HI IZj- They revealed potential gain of such sources over 
systems relying on strongly attenuated laser pulses. However, in these experiments, 
the actual performance of QKD is intrinsically linked to photon statistics of the single 
photon source jH]. 

Following the proposal of De Martini et al. [HI EH], we recently realized a SPS 
based upon pulsed excitation of a single molecule [IT]. Among various experimental 
realizations of single photon sources |T2l CHI UH UE1 UEJ HH UE], a molecular-based 
SPS presents several advantages. First, it can be driven at room temperature with a 
relatively simple setup which achieves global efficiency exceeding 5 % for single photon 
production and detection. Secondly, since the molecular fluorescence lifetime is a few 
nanoseconds, high repetition rate can potentially be used. Finally, background-emitted 
photon intensity level is extremely low, for carefully prepared samples. 

At the single pulse timescale, the figure of merit of a SPS can be characterized by 
efficiency of delivering triggered photons to target and by the ratio of single photon to 
multiphoton pulses [10]. In reference (TTj, we extended this analysis to measurement of 
SPS noise properties over a wide integration timescale range. In the detection scheme, 
complete statistical information is extracted from the "photocount by photocount" 
record. We showed that measured photon statistics strongly deviates from Poisson 
law, therefore clearly exhibiting non-classical features. 

In this article, we detail the steps of this work, from realization of a molecular-based 
SPS to extensive statistical analysis of detected photons. 

2. Single photon emission from a single molecule 

2.1. Principle of the experiment 

As fluorescence light of a 4-level single emitter is antibunched for timescale on the order 
of the excited state radiative lifetime [JjiJ EH EH1 EJJ , such systems can simply produce 
single photons on demand [01113 EH]- As summarized in figure [TJ the molecule is pumped 
into a vibrational excited state by a short excitation pulse. It then quickly decays to 
the first electronic excited state by a non-radiative process [22] • The emission of a single 
fluorescence photon then coincides with radiative de-excitation toward the ground state 
vibrational multiplicity, followed once again by fast non-radiative decay (figure [TJ. To 
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Figure 1. (a) Single- photon generation by pulsed excitation of a single 4-level 
molecular system from the ground singlet state Sq to a vibrationally excited sublevel 
of the singlet state Si. Solid arrows corresponds to optical transitions, whereas 
dashed arrows depict non-radiative fast (ps) de-excitation. In order to emit a single 
photon per excitation pulse, the pulse duration 8t must be much shorter than the 
radiative lifetime 1/r. (b) Absorption and fluorescence spectra of DiICis(3) dye 
embedded in a thin polymer film, measured respectively with absorption spectrometer 
and spectrofluorimeter with 514 nm excitation wavelength. Note that in the SPS 
experiment the 532 nm excitation wavelength is well separated from the dye's 
fluorescence emission which is centered at a wavelength around 570 nm FWHM. 



emit more than one photon at a time, a molecule has to undergo a full excitation, 
emission and reexcitation cycle within the same excitation pulse. The probability of 
this occurrence is extremely small when the pulse duration is much shorter than the 
excited state lifetime [HI Following the theoretical analysis of reference jES, we 
chose a pulse duration of ~ 150 fs, which makes this probability less than 5 x 10~ 5 . 
This value is negligible in comparison to the one associated to parasitic light, such as 
residual fluorescence from the molecular host matrix. To get one fluorescent photon per 
excitation pulse, the repetition period must also be much longer than the excited state 
lifetime so as to ensure relaxation into the ground state before application of the next 
excitation pulse. 

2.2. Experimental setup 

We use standard confocal microscopy techniques to perform selective excitation and 
detection of single- molecule light emission at room temperature [21] . This setup allows 
one to readily achieve two required features for observation of non-classical photon 
statistics, namely good collection efficiency of emitted photons and high rejection of 
optical background noise. 

The laser source, used for fluorescence excitation, is a femtosecond tunable 
Titanium: Sapphire (Ti:Sa) laser, frequency doubled by single-pass propagation into a 
LHO3 nonlinear crystal. The initial repetition rate of 82 MHz is divided by a pulse-picker 
with frequency set to 2.05 MHz (pulsed excitation repetition period r rep = 488 ns) to 
avoid surpassing the maximum electronics counting rate. 

The excitation light, centered at 532 nm, is reflected by a dichroic mirror into an 
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Figure 2. Experimental setup for characterization of the triggered SPS. Single 
molecules are excited by a frequency-doubled femtosecond Ti:Sa laser at 532 nm. 
The laser is followed by PP: pulse picker; C: LHO3 nonlinear \^ crystal; EO: 
ADP(NILtP04) electro-optic cell; P: linear polarizer; PZT: piezoelectric translation 
stage; Obj: oil immersion microscope objective (x60, NA=1.4); DM: dichroic mirror; 
PH: pinhole for confocal detection (30 /im diameter); NF: notch filter centered at 
532 nm; BS: non-polarizing beamsplitter; SPAD: single photon silicium avalanche 
photodiode; TAC: time to amplitude converter; MA: multichannel analyser; TIA: time 
interval analyser (GuideTech, Model GT653) and PC: computer. 



inverted microscope. It is focused on the sample with an oil-immersion high numerical 
aperture objective, leading to a spot diameter of ~ 350 nm FWHM. The fluorescence 
light -redshifted with respect to the excitation- is collected by the same objective, then 
transmitted through the dichroic mirror, and finally focused inside a pinhole for the 
confocal configuration. After recollimation, residual excitation light is removed by an 
holographic notch filter. 

The samples used in our experiment consist of cyanine molecules DiICis(3). This 
dye choice was motivated by its fluorescence efficiency and photostability, with an 
emission spectrum well suited for detection using silicium avalanche photodiode (see 
figure Q (b)). The dye molecules are embedded in a thin layer of PMMA deposited on 
a microscope coverplate by spin coating. The emitters are randomly distributed within 
the PMMA layer (thickness ~ 30 nm) at an approximate concentration of one molecule 
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per 10 fim 2 . 

To ensure localization of a single emitter in the detection volume, we use a 
standard Hanbury Brown and Twiss setup |2S]- It consists of two single-photon- 
counting avalanche photodiodes (SPADs) placed on each side of a 50/50 nonpolarizing 
beamsplitter. A Start-Stop technique with a time-to-amplitude converter allows us to 
build a coincidences histogram as a function of time delay between two consecutive 
photodetections on each side of the beamsplitter. Following the textbook experiment 
of P. Grangier and A. Aspect on quantum properties of single photon states [21], the 
absence of coincidence at zero delay gives clear evidence of single photon emission |12j . 

We hence apply a simple three-steps procedure as explained in reference jTTj. We 
first raster scan the sample at low energy per pulse (~ 0.5 pJ) so as to map the 
fluorescence intensity and locate efficient emitters. We then put the excitation beam 
on a given emitter and measure the autocorrelation function at low excitation energy. 
We hence determine whether a single fluorophore or an aggregate of several molecules 
is excited. Once such preliminary identification has been achieved, the single molecule 
is excited at much higher power so as to ensure saturated emission [TTj . 

2.3. Data acquisition 

Once a single emitter is located, we switch the detection procedure from the Start- 
Stop method to a complete recording of photon arrival times. The properly normalized 
interval function c(r) measured by the Start-Stop technique corresponds to the intensity 
autocorrelation function g^ 2 \r) in the limit of short timescales and low detection 
efficiency |2HJ - However, to characterize more completely the statistical properties of 
the photon stream, one needs to test for correlations on timescales much longer than 
the excitation repetition period r rep . In that case, the relationship between g^(r) and 
c(t) becomes more complicated flFH, I2()j. Instead of solely inferring g^ir) from c(r) 
measurements, we have chosen to keep trace of the full range of dynamics by recording 
every photodetection time with a Time Interval Analyser (TIA) computer board. From 
this set of photocounts moments (that we call timestamps), detected photons statistics 
can then be directly analyzed over a wide range of timescales. Such a procedure avoids 
any mathematical bias in photon statistics analysis. 

The total number of fluorescence photons that can be produced by a single 
molecule is limited at room temperature by its photostability |27| . Under weak CW 
excitation, a molecule of DiICi 8 (3) typically undergoes 10 6 excitation cycles before 
irreversible photobleaching occurs [2H] . In our experiment, the excitation pulses energy 
is progressively ramped up to a maximum value of 5.6 pJ that ensures saturation of the 
So — > Si transition [TT]. This energy ramp, realized using an electro-optic modulator, 
consists of a 50 ms linear rise followed by a plateau lasting 300 ms and linear decrease 
(figure EJ). We experimentally found that applying such a procedure substantially 
improves molecular photostability, compared to an abrupt excitation. We then select 
the timestamps of photocounts that occured during the plateau of the excitation. These 
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Figure 3. Black solid line (right scale): laser pulse energy E p vs. time during 
single molecule excitation. Maximum laser pulse energy i?™ ax of 5.6 pJ saturates 
the molecular transition. Gray solid line (left scale): number of fluorescent photons 
detected during 50 (is integration duration. Photobleaching of the dye occurs 162 ms 
after excitation at the maximum energy per pulse begins, as delimited by two vertical 
dashed lines. 



events correspond to saturated emission of our molecular SPS. Our analysis "photocount 
by photocount" then relies on determination of the number of detected photons in gated 
windows synchronized on Ti:Sa excitation pulses. 

Data is first pre-processed over a discrete time grid. The excitation time base 
is reconstructed from the timestamps ensemble {U} and by applying a time filtering 
procedure described in Appendix A The gate duration is 30 ns, more than ten times the 
typical radiative lifetime of the molecule in PMMA layer. All records outside the time 
gates are rejected, slightly improving the signal to background ratio. Each timestamp ti 
is then attributed to a pulse pi in the time grid, and to each excitation pulse p, a number 
nf^ = 0,1,2 of detected photons is finally associated. The probability distribution of 
the number of detected photons per pulse is deduced from {rip^K as summarized in 
tabled 

We next analyse the photon statistics of a data set extracted from SPS emission 
displayed in figure 01 We have selected photocounts recorded between the plateau 
beginning and molecular photobleaching clearly identified by a sudden drop in 
fluorescence emission. During this time, the molecule was excited at constant maximum 
pumping energy, yielding 15332 photodetection events for 325313 excitation pulses. The 
time filtering mentioned above is then applied keeping 15138 synchronous photocounts 
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Table 1. Single-pulse statistics of a molecular single photon source, as obtained after 
numerical synchronization (see |Appendix A\ . This data will be referred to as (S). The 
total number of excitation pulses in the sequence is 325313, leading to a total of 0,1 or 
2 photon number events of 15138. The mean number of photons detected per pulse is 
(n) = 0.04653. No events with n ( - d - ) > 2 are observed due to deadtime in each detection 
channel. 
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3. Single pulse photon statistics 

The single-pulse statistics presented in table [T] are the direct outcome of photocounts 
acquisition. They correspond to the molecular emission displayed in figure El While 
our SPS photon statistics appear to differ from a classical Poissonian distribution, 
the influence of our experimental setup on these measurements must be considered 
for accurate interpretation of these figures and for comparison with Poisson shotnoise 
reference. 

3.1. Influence of deadtime 

In the following, we make a distinction between the distribution of the photons produced 
by the source, denoted by script notation (V), and photocount statistics, for which we 
conserve the usual notation (P). 

Due to existence of a ~ 280-ns deadtime for each detection channel, a nonlinear 
relationship exists between detected photon statistics and source photon statistics. 
Indeed, for a given excitation pulse, the number of detected photons in a 30-ns gated 
time window cannot exceed two if we use two avalanche photodiodes (APDs) operating 
in the photon counting regime compulsory for our experiment. 

Denoting by V m (n) the photon number probability distribution of incoming light 
on the detection setup, the nonlinear tranformation relating this probability to the 
detected photon probability P{n = 0,1,2) is simply computed for "ideal APDs". By 
"ideal APD" , we mean that each photodiode clicks with 100 % efficiency immediately 
upon receiving a photon, but that no more than one click can occur in a given repetition 
period. In the approach developed here, we consider the ideal APD case for the following 
reasons: 

• limited quantum efficiency of the APD (65 % in our experiment) is included in an 
overall linear loss coefficient along with other linear losses of the detection chain 

• deadtime is shorter than repetition period r rep and much longer than pulse duration. 
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For detection with a single ideal APD, the relationship between the photocount and 
incoming light statistics is 

oo 

P(0) = P in (0) and P(l) = ^"W C 1 ) 

n>l 

With our experimental detection scheme, random splitting of photons on two sides 
of 50/50 beamsplitter gives 

P(0) = P in (0) (2) 

oo 

P(l) = £ P'» ^ (3) 

n>l 

P(2) = £ P in (n) (l - JL) (4) 

n>2 ^ ' 

The relationship between Pin) and photon statistics P(n) in SPS emission, comes 
from accounting for linear attenuation between SPS and detection. We call r\ the overall 
detection efficiency, which includes all linear propagation losses and photodetector 
quantum efficiency. V m is then related to V by the following binomial law 

oo 

r' m (n) = £ O (1 - rjy™ V(m) (5) 

m=n 

Combination of equations ©, ©, © and © leads to a direct analytical relation 
between V(n) and P(n = 0, 1, 2). 

Note that existence of such a saturation limit, due to detection deadtime, has no 
influence on photocount statistics of a perfect SPS, for which V{n > 2) = 0, as long as 
excitation repetition period is longer than electronics deadtime. On the contrary, for 
a "real" source with background light, the number of detected multi-photon pulses is 
systematically underestimated, leading to statistics artificially squeezed in comparison 
to shotnoise reference. 



3.2. Calibration with a coherent source 

3.2.1. Coherent beam photocount statistics SPS performance can be directly evaluated 
by comparing single pulse photon statistics with those of a coherent source. This 
calibration takes into account the linear and nonlinear effects of our detection setup 
and permits accurate measurements of the multi-photon events probability reduction 
between single photon and Poissonian sources. 

The photon number probability distribution for a coherent pulsed beam (C) is given 
by a Poisson law. According to equation ©, linear loss between the source and APDs 
change the mean photon number per pulse a to rja while the photon statistics remain 
Poissonian. The expected photocount statistics can then be calculated by applying 
nonlinear transformations (equations (J2J), © and to a Poissonian distribution of 
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Table 2. Photocount probabilities P{n) for SPS (S), reference experimental coherent 
source (R), and theoretical coherent source (C), for which photocount statistics are 
affected by detection. This table also displays the mean number (n) of detected photons 
per pulse. 





n = 1 


n = 2 


(n) 


Ps(n) 


0.04644 


4.6 x 10" 5 


0.04653 


P R (n) 


0.04520 


50 x 10" 5 


0.04620 


Pc(n) 


0.04514 


53 x 10" 5 


0.04620 



parameter rja 

P c (0) = e-« a (6) 

P c (l) = 2e-" a/2 (l - e-" a/2 ) (7) 

P c (2) = (1 - e-^ 2 )\ (8) 

such distribution being termed as Pc(n) in tabled 

3.2.2. Experimental calibration A strongly attenuated pulsed laser beam is used as 
experimental reference to mimic a pulsed coherent source. It is obtained by slightly 
detuning the Ti:Sa wavelength from the notch filter rejection band resulting in detection 
of residual pump light reflected from the sample. The photocount statistics of this 
experimental reference (R) are then compared both with the experimental single photon 
source (S) and the calculated photocount distribution expected from a Poissonian source 
(C). To establish a valid comparison, calculated and experimental calibrations are 
determined for an -almost- identical mean number of photons detected per pulse. 

Table El shows that theoretical predictions for the coherent source are in good 
agreement with experimental calibrations, proving that our detection model accounts 
for all significant biases. We can therefore confidently interpret the molecular SPS 
photon statistics we measure. For our SPS, the number of two-photons pulses is 10 times 
smaller than the corresponding probability for a Poissonian source. As mentioned earlier, 
residual multi-photon pulses mostly results from background fluorescence light triggered 
by Ti:Sa excitation. Indeed, the wavelength of this parasitic light lies within the 
molecule's fluorescence band and therefore cannot be filtered out. Careful optimization 
of substrate purity as well as of the chemicals purity used in sample fabrication can 
likely lower background fluorescence. 

3.3. Molecular SPS efficiency 

Our molecular SPS emission can be modeled as the superposition of a perfect 
SPS and a coherent state of light. In this model, all sources of linear loss 
(production + collection + detection), are gathered as an overall efficiency rj. For 
this perfect SPS (perf. SPS), the photon probability distribution of light inpinging on 
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the APD is then given by 

*£rf.sra(0) =1 ~V 

^erf.SP S (l) =V (9) 

Ka.spsin > 2) = 0. 

Background (backgnd.) emission is modeled by a coherent state of light with a mean 
number 777 of detected photons per pulse. The corresponding photon probability 
distribution is then 

TTackgnd. (n) = ;, /7J , for n > 0. (10) 

Applying equations Q to (JH) to the (perf. SPS + backgnd.) probability 
distribution leads to the following analytical expressions for the real single photon source 
(S) photocounts statistics : 

P s (0) = e-™(l-rj) 

P s (l) = 2 {e~ m/2 - e~ rn ) + n {2e' rn - e~ rn/2 ) (11) 
P s (2) = (1 - e -^ /2 f + 77 (e"^ 2 - e"" 7 ). 

Values for collection efficiency 77 and signal-to-background ratio I/7 can be inferred 
from measured photocount statistics Ps (see table EJ. Using equations (JTTJ) for 
experimental values of Ps(l) and Ps(2), we finds 77 ~ 0.04456 and 777 ~ 2.02 x 10~ 3 . 
This leads to a signal-to-background ratio of 22, in good agreement with that measured 
by sample raster scan. 



3.4- Single-pulse Mandel parameter 

From a statistical point of view there exists two main differences between experimental 
and ideal SPS: source overall efficiency lower than unity, and finite ratio of single-photon 
to multi-photon pulses. Light produced by an ideal SPS consists in the periodic emission 
and detection of single photons with 100% efficiency, its intensity fluctuations being then 
perfectly squeezed. On the other hand, a real SPS yields less squeezing [Bj. 

It is then meaningful to assess SPS performance by measuring its intensity noise 
on the excitation repetition period r rep timescale [TT]. Such analysis requires evaluation 
of the single pulse Mandel parameter Q [29J. This parameter characterizes deviation 
of photon statistics from Poissonian statistics for which Q = 0. Subpoissonian (resp. 
superpoissonian) statistics correspond to Q < (resp. Q > 0). For the distribution 
{rip d ' ) } of detected photon number, the Mandel parameter is defined by 

(n) (n) 

where (n) stands for the average value of {fip^} calculated over the ensemble {p} 
of excitation pulses. Note that an ideal SPS would yield Q — —1. Moreover, for 
any statistical distribution, the effect of linear attenuation can be straightforwardly 
evaluated: after linear attenuation 77, a Mandel parameter Q would be changed in r]Q . 
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This means that every statistical distribution converges towards Poissonian statistics 
under attenuation. This sensitivity to loss for measurements of non-zero Mandel 
parameters is similar to sensitivity observed in squeezing experiments that measure 
reduced photocurrent noise spectra with respect to shotnoise reference. 

For our molecular SPS, the Mandel parameter Q of the photocount statistics can 
be computed directly from single-pulse photocount probabilities 



From table 121 data we infer a Mandel parameter Q$ = —0.04455 for the SPS. This 
negative value for Q confirms that our SPS indeed exhibits subpoissonnian statistics at 
the timescale r rep . Since very few multi-photon events are observed, the value of Qs 
is almost only limited by the collection efficiency which imposes a lower limit on Q: 



Our measurement of Qs can then be compared to a Poissonian reference 
measurement. Here again, statistical bias introduced by APD deadtime must be 
taken into account. From equations ©-© and (fT3j) . we can derive the Mandel 
parameter of detected photons for a coherent source (C) of parameter a. Noticing 
that (n) c = 2(1- e~ Q/2 ), we have 



As a consequence of photodetector deadtimes, a coherent source gives subpoissonian 
distribution of photodetection events. In our case, a coherent source with the same mean 
number (n)c = (n) = 0.04653 of detected photons per pulse as the SPS, would then 
yield Qc = —0.02327 > Qs- Despite this detection bias, our direct measurement of the 
Mandel parameter, still yields a value for Qs that clearly departs from that of Poissonian 
statistics. This measured Mandel parameter is larger (in absolute value) than those 
measured in previous measurements by more than one order of magnitude [2HJ 12 120] • 

4. Single photon source intensity fluctuations 

Emission intermittency has been observed with most single photon sources realized 
so far [ini E3 HE! • This effect decreases source efficiency and contributes to additional 
source of noise. Better understanding of physical processes responsible for intermittency 
would likely lead to significant improvement of current SPS devices. 

To characterize intermittency for our molecular SPS, we have investigated its 
influence on the photon statistics recorded with the time-resolved photon counting 
system. For a periodically trigerred SPS, this analysis is equivalent to study of source 
intensity noise over a wide range of timescales, which is usually done in the frequency 
domain for squeezing experiments, using a radio-frequency spectrum analyzer. 




(13) 



Qnmit = -V = -0.04456. 




(14) 
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4-1. Measuring intensity fluctuations : Time-varying Mandel parameter Q(T) 

The analysis performed on single-pulse photon statistics (see section 13. 4 J) can be 
extended to multiple-pulses scale, allowing characterization of intensity fluctuations 
at any timescale greater than the pulse repetition period. To do this, we analyze 
fluctuations of the total number N (T) of photons detected during an integration time 
T = A^-r rep , which is a multiple of the repetition period. This analysis therefore 
corresponds to study of statistics of the photocounts number recorded during M. 
successive pulses. 

We then introduce the time- varying Mandel parameter Q(T) (20] , defined similarly 
to single-pulse Mandel parameter. To perform statistical analysis, we extend the 
procedure used in section El More precisely, the complete data {rip } p= i y ,^ 
corresponding to photocounts recorded during M consecutive excitation pulses is split 
in successive samples, each lasting T. We then obtain ACampie = E [M /M] samples. We 
call Nk (T) the number of photocounts recorded during the k th sample. We then have 

r {k+l)T (k+l)M-l 

N k (T)= / I(t)dt= J2 u p- ( 15 ) 

JkT P =kM 

The statistical average over these samples of duration T is denoted ( )t, and we hence 
have 

-A/samplc 1 

Jv sample fc=Q 

Using this notations, the time-dependent Mandel parameter is given by 

that allows direct comparison of SPS noise properties to those of Poissonian ligth beam. 



4-2. Intensity noise and intermittency in the molecular fluorescence: the ON-OFF 
model 

To analyze our experimental results and link them to physical parameters of molecular 
fluorescence, we use a simple analytical model of molecular intermittency, in which we 
assume that the SPS can be ON or OFF. We call p the ON to OFF transition rate and 
q the OFF to ON one. These rates correspond to lifetimes r on = 1/p and r off — 1/q, 
respectively. 

4-2.1. Physical interpretation of the ON-OFF model For molecular SPS ^Tj and other 
SPS's relying on fluorescence of a single emitter (e.g. single NV centers in diamond 
nanocrystals JH]), ON-OFF intermittency stems from the presence of a metastable non- 
fluorescent excited state in the energy level structure. Dynamics of ON-OFF behavior 
can then be computed from the three-level structure shown in figure |U 
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Figure 4. Three states energy level structure and the corresponding ON and OFF 
states in the SPS intermittency model. In the ON state, the molecule undergoes 
fluorescent cycles between the ground So and excited S\ singlet states. In the OFF 
state, the molecule is trapped in the dark metastable triplet T\ state. Coupling from 
S\ to T\ occurs at each excitation pulse with the intersystem crossing probability PisC) 
yielding a transition ON— >OFF rate p — 'Pisc/ T rcp- The reverse transition OFF— >ON 
occurs at rate q — 1/tt, where tt is the triplet state lifetime. 



• ON — > OFF transition corresponds to relaxation from the optical excited state 
Si to the triplet state T\. For each excitation cycle, the probability Vise of this 
intersystem crossing process is very small in the case of DiICi8(3) molecule used 
in our experiment {Vise — 10~ 4 ). Moreover, since singlet-triplet transitions occur 
exclusively from the excited state Si, the excitation repetition period must be 
considered in defining the source ON state lifetime T on , which is then r on = T rep /Visc, 
assuming saturated excitation regime. 

• OFF — ► ON transition consists simply of non-radiative decay from triplet T\ to 
ground Sq state. Note that the triplet level is metastable since selections rules 
forbid direct optical transition to the ground state. The triplet state lifetime 
tt = T Q fj = l/q is therefore usually much longer than a typical fluorescent lifetime 
(in the case of DiICi 8 (3), r T ~ 200 y& |2HJ ) - 

4-2.2. Dynamics of the ON-OFF system Under periodic pulsed excitation, the ON- 
OFF dynamics can be described using a discrete time model. As transitions between 
ON and OFF states are random, we introduce a stochastic variable to account for 
the source state at instants tk = kr Tcp . This parameter has value = 1 (resp. = 0) 
if the source is in the ON (resp. OFF) state at time tk- 

We then call the probability for the source to be in the ON state at time t^. As 
the SPS emits photons exclusively from the ON state and never from the OFF state, 
also corresponds to the photoemission probability at time We assume that lifetimes 
r on and T g of the ON-OFF states are much larger than the repetition rate r rcp . Then, 
the ON— >OFF transition probability is pr rep and the OFF— >ON transition is qr rcp . It 
follows that the state of the emitter at pulse k + 1 depends only on its state at pulse k. 
The recursion relation for the probability Uk+i of the source to be ON at time tk+i is 

u k+ i = (1 - pr Tcp ) u k + qr Tcp (1 - u k ) , (18) 
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which leads to the general solution 

u k =(u — J (1 -pr rcp - gr rcp ) fc H — . (19) 

V P + qj P + q 

Stationary probabilities for the molecule to be either ON or OFF are then 

P - = ¥TT) (20) 
p °" = l - p °« = vky (21) 

4-2.3. Source intensity and Mandel parameter vs. time According to our model, light 
emitted by the source is a succession of single photon pulses emitted at time tk = kr Tep 
with probability Uk, corresponding to intensity 

+ 0O 

I(t) = 2J <^(* ~ kr rcp ) x r k , with r k = or 1. (22) 

k = —00 

The recursive relation (|T9|) for the ON-OFF model permits computation of 
statistical properties of source intensity I(t). In particular, we can derive the time- 
dependent Mandel parameter from the variance of the number N(T) of photons emitted 
by the intermittent source during T = A4.T vep . Details of this calculation are given in 
|Appendix B| 

Analytical expression of the Mandel parameter given by equation (|B.13jl . can be 
simplified in the regime for which (5 = (p + q)T Tep <C 1, leading to the following Mandel 
parameter expression for a "perfect" SPS with intermittency 

g pcr ,sPs(-Mr rcp ) = 2 -Pjp^ {l - JL [1 _ (i _ ft"] } - I. (23) 

Experimental measurements of the Mandel parameter are also affected by overall 
efficiency 77 smaller than unity (see section I33|) . Taking into account this limitation 
which is equivalent to linear loss, the Mandel parameter of the real source Qs(T) is 
given by 

Q S (T) = r / g perf . S Ps(T). (24) 



4.2.4. Experimental data analysis As shown in figure EJ our experimental data are 
well-fitted by equations (f2*3"|) and over more than four orders of magnitude in time. 
Setting the measured efficiency to r\ — 0.04456, the fit yields pr rep = Vise — 2.1 x 10~ 4 
and tt = 250 fis, for the remaining two free parameters. These values are in good 
agreement with values given in reference [28J. 

Figure clearly shows that source photon statistics differ on short and long 
timescales. On timescales shorter than ~ 8r rep , the Mandel parameter of the source 
Q(T) is smaller than Qc, the theoretical value of the Mandel parameter for poissonian 
light including the detection deadtime (horizontal dashed line on figJHJ). On this short 
time scale, the SPS's photocount statistics are those of non-classical light. On timescales 
larger than ~ lOfis, fluorescence intermittency due to the triplet state, influences the 
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1 0" 6 1 0" 5 

Integration time 7~(s) 

Figure 5. Direct measurement of Mandel parameter Q(T) over short integration 
time T. The dashed horizontal line shows Q{T) for the equivalent coherent source 
(C), taking into account detection dead time. Inset shows Q(T) for longer integration 
time. The solid curve is a fit given by the model accounting for intermittency in SPS 
emission. 

photocount statistics by introducing excess of noise resulting in a positive value of the 
Mandel parameter. 

The model developed here for a perfect intermittent SPS fits our experimental 
data with good accuracy. Indeed, apart from detection loss, other imperfections can be 
ignored or handle by the following: 

• since the repetition period r rep is much longer than the photodetection deadtime 
and since multi-photon events are extremely rare with our SPS, APD deadtime does 
not alter significantly the photocounts statistics. The detection can be considered 
effectively linear and equation remains valid in the presence of detection 
deadtimes. 

• high signal-to-background ratio means that background light does not contribute 
significantly to photocount statistics. It can therefore be neglected, as done 
implicitly in the model developped in this section. It can moreover be shown 
that addition of uncorrelated Poissonian background light of intensity B to the 
perfect SPS signal S is equivalent to loss. If we model the real source by the 
superposition of fluorescence background and light from a perfect SPS, then, 
introducing p = S/(S + B), the Mandel parameter of the real source is simply 
given by Q s+B = pQ S - 

4-2.5. SPS intensity autocorrelation function As a consistency check for our study, the 
time dependent Mandel parameter analysis can be compared to a different approach 
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Time delay (x x rep ) 



Figure 6. Discrete-time photocount autocorrelation function , computed from the 
set of data {rip } studied, in the range of time delay T = 1 x r rop to T — 1000 x r rcp . 
The dashed horizontal line corresponds to shotnoise reference value. 



using the intensity autocorrelation function pH] . the measurement of which being 
at the heart of fluorescence correlation spectroscopy 

In terms of discrete time variables, a discrete time autocorrelation function for time 
delay A x r rop is given by 



where n» is the number of detected photons in the i excitation pulse, and A is an 
integer. 

This discrete correlation function is directly related to the intensity autocorrelation 
function g^(r) = (I(t)I(t + r))/(/(t)) 2 usually measured with Start-Stop 
techniques |12j . It can be shown that the normalized area of the k th peak (with the 
gth re f erence peak corresponding to r = 0) of g^(r) over a period r rep is equal to 



The ON-OFF model developed in section |U can be applied to calculate G^(A) 



Note that the latter value is chosen because blinking occurs in a timescale range of 
~ 1000 x r rep . Results of this G^(A) calculation are displayed in figure El The 
experimental curve is fitted with equation providing another way of measuring 

dynamical parameters of intermittent molecular SPS. This fit yields Pise — 1-6 x 10 -4 
and tt = 180 fis, which are in good agreements with the values obtained in section 14.2.41 
using Mandel parameter analysis. 

Note that, on short time scale, the statistical noise is higher on than on Q. 
This is due to the fact that is computed over fewer but bigger statistical samples. 



G^(A) 




(25) 



G( 2 )(A = k). 
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5. Conclusion 

We have realized an efficient triggered single photon source relying on the temporal 
control of a single molecule fluorescence. After a comparison to Poissoninan coherent 
light pulses with the same mean number of photons per pulse, we have characterized 
intensity noise properties of this SPS in the time domain and photocounting regime. 

From the record of every photocount timestamp, we calculate the second order 
correlation function or equivalently the time- dependant Mandel parameter Q(T). 
Observed negative Q(T) values signifie non-classical photocount statistics. 

This time-domain analysis is complementary to fluorescent correlation spectrocopy 
techniques for investigating photochemical properties at the single-emitter level. More 
specifically, we have modeled fluorescence intermittency by a two-state ON<->OFF 
dynamical process. By fitting a theoretical analytical expression of the Mandel 
parameter for an intermittent SPS, we obtained quantitative values for relevant 
molecular photodynamical parameters. Such a direct time-domain statistical analysis 
could give insight into molecular properties such as conformational changes [33J • 
resonant energy transfer or collective emission effects in multichromophoric 

systems |3E]- 

With expected application to quantum cryptography, higher overall efficiency 
within a given emission spectral band should be reached so that single photon sources can 
exhibit advantages over attenuated laser pulses [Hj. In recent experiments we coupled the 
fluorescence of a single emitter (a colored center in a diamond nanocrystal) to the single 
mode of a planar microcavity and observed a significant increase in spectral density of 
the emitted photons. These preliminary results are promising realization of an efficient 
single photon source well-suited for open-air quantum key distribution. 
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Appendix A. General analysis technique of a set of photocounts 

A set of data consists of a list of timestamps recorded by the Time Interval Analyser 
computer board. In this appendix we describe the protocol developed to process raw 
data. This procedure allows us first to postsynchronize the timestamps on an excitation 
timebase and then to build the set {rip } of the number of detected photons for each 
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excitation pulse p. 

The pulsed excitation laser acts as a periodic trigger of emitted photons with 
repetition period r rcp ~ 488 ns. An excitation laser pulse is emitted at time t s tart+pxr rep , 
where the pulse is indexed by the integer p. The parameter t start represents the pulse 
emission time taken as the first (p = 0) of the data. For each set of data, t s tart arid r rep 
must be determined because the repetition period of the Ti:Sa femtosecond laser can 
fluctuate slightly between acquisitions. However, the laser repetition rate is stable over 
the typical acquisition duration (under one second), and r rep is therefore constant for a 
given data record. 

Single photon emission by the molecule occurs at each excitation pulse after a 
random time delay related to the molecule's excited state lifetime. Non-synchronous 
photocounts due to APD dark counts are rare, so almost all recorded photocounts 
are triggered by photons emitted by the molecule, with few by photons from residual 
fluorescence background. For these reasons, t sta rt and r rep can then be determined 
directly from recorded data. 

The i th photocount timestamp ti can be expressed as 

U = tstart + (pi x r rep ) + 5n, (A.l) 

where integer pi e {l,...,Af} indexes the laser pulse preceding detection of the i th 
photon, and data to be analysed lasts N repetition periods; Sr^ is the time delay 
between excitation pulse and photocount timestamp (0 < 5ti < r rop ). Given a set 
of timestamps {ti}, relevant information can be equivalently represented by lists {p^ 
and {S^}, provided the value of r rep is known accurately enough. Since the fluorescence 
lifetime of the molecule is much shorter than the laser repetition period, Sri <C r rep , as 
long as photocount i is not a dark count, which is rarely the case. 

As the laser period r rep is not known precisely, we first attempt to synchronize the 
data on the excitation timebase considered as a clock of period T c i ock close to the expected 
laser period. We introduce a delay function parametrized with i start and r c i oc k, that gives 
for each timestamp ti the time delay between this timestamp and the corresponding top 
of the clock + 

Dela yi st a rt ,r clock (^) =u- Wt - e ( \ d ^ art ) x Tclock ' ^ A ' 2 ) 

where, if the clock period differs from the laser period, the drift of the time delay baseline 
with the pulse index pi is linear as ican be seen in figure lAlf a) 

Dela yt st a rt ,r clock (^) - Pi( T rep ~ T c lock) - ~^ Orep ~ T clock ) . (A.3) 

'rep 

From the slope of the delay function baseline, we infer a new value for r c i oc k. Note that 
the first guess for r c i oc k is usually so far from the laser period that the delay function 
takes a saw-toothed shape, each jump corresponding to the delay reaching a multiple 
value of r c i oc k. As a consequence, only a linear fraction of the sample corresponding to 

+ EQ stands for the integer part function. 
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Figure Al. Synchronization procedure of timestamps {ii}, showing delay function 
Delay(ti) vs. pulse index E[{ti — i s tart)/ T ciock] for a given set of parameters i s tart an d 
Tciock- (a) Case with a small linear drift of delay baseline, when r c i oc k is close to T rcp 



but tstart is incorrect, (b) Case when T c i oc k 
istart is properly estimated. 



T rep within relative precision of 10 and 



a single saw tooth can be used at first. In further steps, estimation of r rep improves and 
fewer jumps occur. Longer samples, corresponding to higher fit precision can then be 
processed. This procedure is repeated until the whole data set is used, leading to the 
situation of figure IXTf b). It corresponds to the same fraction of data as in figure IXlT a). 
for which r rcp is determined up to relative precision greater than 10~ 9 . 

Once T rep and t sta rt values are known, calculation of the lists {p^} and {6^} is 
straightforward using 

Pi = E (t-i^^) and 5r t = Delay W)Tiep (^). (A.4) 

A time filtering procedure is then used to eliminate all photocounts with time delay 
much longer than the molecule excited state lifetime. To implement this filter, we use 
a time window of duration AT window . From the set {pi}, we calculate the number of 
photons detected by the two photodiodes in the time interval [piT rep , PiT rep + AT win d ow ]- 
The time window duration AT wilK i ow must be shorter than the laser period and 
much longer than the molecular excited state lifetime 1/T so that the probability of 
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discarding a "real" photodetection event is negligible. The chosen time window duration 
AT window = 30 ns meets these two conditions, considering 1/T ~ 2.5 ns for DiICi 8 (3) 
dye. Note that choice of a time window significantly shorter than the laser period has the 
advantage of filtering out our data from the majority of non- synchronous background 
photocounts, such as APD dark counts. 

The processed data, now expressed as the table {rij,pj}, and shortened to {rip , } in 
the body of this article, allows us to characterize the statistics of our source on timescales 
from T rcD ~ 500 ns to milliseconds. 



Appendix B. Statistical characterisation of an intermittent SPS 

In this Appendix, we derive a general analytical expression of the "perfect" intermittent 
SPS Mandel parameter Q(T) defined by equation (fTTj) using the ON-OFF model 
introduced in section [Ol We also retrieve the approximate expression (|2^j) for Q(T). 
We assume that source emission has reached its steady state at time t = . 

The total number N (T) of photocounts recorded during an integration time 
T = M -Trcp corresponding to M. consecutive excitation pulses is related to the stochastic 
photocount variable r k associated to k th excitation pulse (see section 14.2. 3|) by 

M-l 

A(T)=]Tr fc , (B.l) 

k=0 

Calulating Q(T) is equivalent to evaluating the variance {N 2 ) T - (N) T 2 , where mean 
values ( )t are defined by equation (fTp^) . Both (N) T and (N 2 ) T should then be evaluated. 
The mean value of N(T) is given by: 

M-l 

<A0 T =£(r fc >=.M-^-, (B.2) 

where the steady state expression of (r k ) = q/(p + q) comes from the definition of 
and the recursive law (|T^|) . Similarly, (N 2 ) T follows 

M-l M-l 
k=0 k'=0 

Index change I = \k — k'\ in the previous equation then yields 

M-l M-lM-k 

k=o k=o e=i 

M-l M-l 

= ^(r fc )+2^(A<-£)C(£), (B.4) 

k=0 1=1 

where we have introduced the discrete time correlation junction C (£) defined as 

C(£) = (r k r k+e )-(r k )(r k+i ). (B.5) 

Calculation of Q(T) now relies on evaluation of C(£). We recall that, in the model 
of section |4~21 source dynamic is described by stochastic process u k , the probability for 
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the molecular system to be in the ON state. The general expression of u k+ e follows from 
the recursive law (fTTJj) 

Uk ; — (1 -pr Icp - gr rep ) H ■ — (B.6) 

p + qj p + q 

The stochastic variable r k equals 1 if the state is ON with probability u k , and if 
the state is OFF with probability 1 — u k - The product r k r k+ g is then equal to 1 only if 
both r k and r k+ £ are simultaneously equal to unity and otherwise equal to zero. It can 
be summarized as 

(r k r k+e ) = P(r k = l)P{r k+t = l\r k = 1) (B.7) 
(r k ) (r h+t ) = P{r k = l)P{r k+i = 1). (B.8) 

Note that P(r k = 1) is the probability that r k = 1, and P(r k+ e = l\r k = 1) is the 
conditional probability for r k+ £ = 1 when r k — 1. To fulfill this later condition r k = 1, 
one needs to have u k = 1. Moreover, by definition, the steady state probability is 
P{r k = 1) — q/(p + q), so that we have 

P{r k+e = l\r k = 1) = -^-(1 - pr vcp - qr rcp f + (B.9) 

p + q p + q 

and, as a consequence, 

c ® = (^T^ (1 " pTrcp ~ qTvep)e - (B - 10) 

This value C(£) is then introduced in equation (|B.4|) . and the expression for the variance 
follows from (EH) and (ED. 



(N 2 )t - (N)> ~ Pq 



M 2a 



M 



1 — a ( 1 — a) 



(B.ll) 



(P + 

where at= (1 - p r rcp - q r rcp ) . 

The general analytical expression of the "perfect" intermittent SPS Mandel 
parameter is finally deduced by 

(N 2 w - m 

{N) T 

p (2- (3 2(1-/3) 1-(1-/3) M ' 



g p erf.sps(-Mr rep ) = V ^— V /r - 1 (B.12) 



P + q\ p m p y" 1 ' (R13) 

where (3 = (p + q) r rop . In the limit of /3 <C 1, which is the case for molecular system 
dynamics considered in the body of this article, we retrieve expression 
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